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A Note on the Stability of Cable-in-Conduit Conductors with Current Dependent 

Power-Law Index 
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A new stability criterion for cable-in-conduit superconductors with power-law current-voltage 
characteristic and current dependent power-law index, n = n{Ic) is given. After a short discussion 
of the power-law volt-ampere characteristic, different models of stability are discussed with spe- 
cial stress on the differences to the older stability models. A typical stability case, the extended 
cryostability is analyzed in detail. This model is characterized by a smooth superconducting to 
normal transition and a power-law type heat generation. The change in the helium temperature, 
typical for a cable-in-conduit conductor, is included in the calculation. Finally, the connection and 
interrelation aspects between the two common approximations of the voltage-current characteristic 
of the technical superconductors, the power-law and the exponential function are investigated. It is 
shown that the exponential form is incompatible with the power-law form if the power-law index is 
a function of temperature and magnetic field. An alternative exponential form is proposed. 
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I. INTRODUCTION 

Because the energy stored in a superconducting mag- 
net, both magnetic and mechanical, can easily be con- 
verted into heat, upsetting the thermal equilibrium of the 
winding, the superconductor temperature is difficult to 
be controlled and the stability issue of superconductors 
became a crucial matter. Therefore, the complete under- 
standing of the thermal behavior of a superconducting 
cable is of great importance both theoretically and ex- 
perimentally. 

A crucial role in the stability issue is played by the 
way the heat generation in the conductor takes place i.e. 
by the volt- ampere characteristic (VAC) of the super- 
conducting material. With the advent of ceramic high 
temperature superconductors and the ITER choice to use 
Cable-in-Conduit Conductors (CICC) for all its magnets, 
the attention on the old non-linear VAC known as Power- 
Law Conductor (PLC) has been refocused. It is expressed 
analytically by the equation 



E ^ Er 



(1) 



where E and / are the electrical field and the operating 
current and Ec is a "man-made" ad-hoc voltage criterion 
used to define the critical current Ic- The critical current 
is itself a function of conductor temperature Tcond and 
magnetic field B. An important role is played here by the 
parameter n, appearing in the exponent. It is responsi- 
ble for the non-linearity of the volt-ampere characteristic 
and the standing heat generation in the real supercon- 
ductors. With this model one can describe almost every- 
thing between a normal conductor {n = 1) and a perfect 
superconductor (n — >■ oo). If n —> oo the volt-ampere 
characteristic of "classical Bean-model" superconductors 
is recovered. For n = 1 the classical resistor (Ohm's law) 



is obtained. 

Based on experimental results we are thought that n 
is by no means always "very high" and it is by no means 
constant. In high temperature superconducting (HTSC) 
tapes it is maximum 15-20 but also values as low as 5 
are frequently reported in the literature. For low tem- 
perature superconductors (LTS) n can be very low (3-5) 
or medium high (15-25) for NbsSn conductors. Slightly 
higher values are reported for NbTi based conductors. It 
is important to distinguish between the strand n-value 
and the cable n- value. There are only few cases when the 
two coincide. The cable n-values are systematically lower 
that the strand values but there are also exceptions to the 
rule. In some NbTi cables with low current transfer ca- 
pacity (insulated strands) the cable 7i-value was higher 
that in the strands used to manufacture the cable [l], Q] . 

The last issue concerning the index of the power law 
which is probably known for a long time but was only 
recently recognized as important, is the dependence of n 
on temperature and magnetic field. This dependence can 
be expressed in almost all cases as a dependence on only 
one variable, the critical current Ic = Ic (Tcond, B) and 
is characteristic for both single strands and cables. For 
NbsSn an additional dependence on strain appears. 

The physical origin of this dependence is still obscure. 
Sure is that n(7c) is large at large Ic (low temperature, 
low field or both) and low at low Ic (higher temperature, 
field or both). Since both limits are well defined ther- 
mally and magnetically it is not possible to assume any 
other hidden factor which could explain this property. 
Concerning the functional dependence of n on /c, lin- 
ear, polynomial and power-law dependencies have been 
reported or assumed. A dependence of the form 



n{Ic)^kir (2) 

with m = 0.3 — 0.5 and a linear dependence 



. (/c) = a + bic 



(3) 



will be considered here. They cover more or less much 
of the spectrum of dependencies found in the literature. 

The paper is organized as follows. In II we first re- 
view different one-dimensional (ID) stability models and 
investigate the impact of changing from the composite 
parallel circuit model of heat generation to the heat gen- 
eration by the power-law (also known as index heating). 
In III a stability model of the unconditional type with 
a PLC-type VAC is investigated. Finally, in IV the re- 
lation between the PLC with variable n-index and the 
exponential form of VAC is investigated. 



II. MODELS OF STABILITY 

With the reference to a superconducting strand, the 
following constituents are relevant for stability. The con- 
ductor itself, composed of superconducting filaments and 
the stabilizing matrix (copper), and the coolant (liquid 
or supercritical He for LTS, LN2 for HTSC). 

If we adopt for simplicity a ID model, the conductor 
Tcond and helium The temperatures are governed by the 
following equations 



J 



P cond^ cond-^ 



dTe. 



d 



iAr 



dT, 



cond 



phe^he-^he 



at ax \ ox 

\^ \-^ cond}^) '^Pw \-^ cond -^he) "T~ ^pulse {J'J 

dTu 



dt 



rhChe—^ = hpw (Tc, 



The 



(4) 



where G {Tcond, t) ^ E {Tcond, t) I {t) is the heat gener- 
ated in the superconductor and Ppuise {t) is the pulse cre- 
ating a deviation from equilibrium. The transient is ini- 
tiated either by the time dependence of the electric field, 
of the current or from an external heat pulse. Depend- 
ing on what terms are retained in the above equations 
and how the electrical field is represented as a function 
of current, different models of stability can be defined. 
This is illustrated in Table 1 and 2 for the conductor and 
helium temperature equations. The terms taken into ac- 
count in each case are marked by the x sign and the 
resulting model name is given in the last column. Let us 
start with the case when we completely neglect the he- 
lium equation. We assume that the helium temperature 
is constant The = Tb and adopt the parallel circuit model 
for the heat generation in the composite [3| i.e. the heat 
generation is formulated as 



G°° {Tcond, t) — < 
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for Tcs < Tcond < Tc 
for Tcond > Tc 



(5) 



In this case we recover the whole class of old models 
of stability for bath cooled magnets and " perfect super- 
conductors" with n — oo. 

For instance, if we assume stationary conditions by 
setting the time derivatives in Table 1 and Table 2 to 
zero we recover the classical " cryostability" or "Stekly" 
model i, [1] 



a 



G°° {Tcond) = hpu, {Tcond ^ Tb) (6) 

the " equal area" model of Maddock, James and Norris 



— ( KAcond 7^ 1 + G^ {Tcond) — hpw {Tcond ~ Tb) 

(7) 



or the "minimum 
Martinelli and Wipf 



ropagating zone" (MPZ) model. 
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by appropriately choosing two of the three terms left 
in Eq.g]). 

Most of the properties of the above models are due to 
a particularity of the functional form of the heat gener- 
ation term as expressed by Eq.([5l). As can be seen, the 
heat generation in the strand is exactly zero below a cer- 
tain temperature, Tcs -the current-sharing temperature. 
With the finite n conductors, where G is non zero at all 
temperatures, this property does not hold and drastically 
changes are expected. A characteristic of power-law con- 
ductor models is that there is always power generation 
at all temperatures. 



TABLE I: Models for stability, conductor equations. 



Heat Accumulation 



Heat Conduction 



Internal Heat Generation 



Heat Exchange 



Model 
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Cryostability 

Equal area 

MPZ 

Dynamic stability 

Protection 

Adiabatic NZP 



TABLE II: Models of stability, helium equations. 



Heat Accumulation 


Convection 


Heat Exchange 


Constant Temperature 


Model 
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X 
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X 
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X 


Pool cooling 

Generalized transient 

CICC 



The cryostability model Eq.®, has an additional par- 
ticularity. Due to the disappearance of the conduction 
term, it is not a differential equation anymore and only 
local values are involved. It can be derived directly from 
a OD model by neglecting the time derivative. 

If the time dependence is not neglected, Eq.(|4l) gives 
other two models for stability. First we have the dynamic 
stability model of Hart ^ 



Pcond^cond^cond rx, 
Ot 



dTcond _ '^ ( A dTcond 

+ G (Tcond, t) ~ hpnj (Tcond ~ The) 

where practically all terms are considered, then the so 
called "protection" model Q 
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which completely neglects cooling and conduction ef- 
fects. Finally we have also the adiabatic normal zone 
propagation (NZP) model [lO| 
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where only the cooling effect is neglected. 

For CICC there is a finite amount of helium at each 
cross section of the cable and the second part of Eq.Q 
cannot be neglected. Also heat generation in CICC is 
better described by power-law. With these too facts in 



mind we can start building stability models for CICC 
picking terms in Table 1 and 2. Again, first assuming 
stationary conditions, we have the following stationary 
models: 

a) Extended Cryostability Model 
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b) Extended Equal Area Model 
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In the class of time dependent models we can imagine 
the following case 

c) OD Transient Stability Model 
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but as discussed above it make sense only for G = G°° . 



III. SMOOTH SUPERCONDUCTING TO 
NORMAL TRANSITION 

In this section we will concentrate mainly on the ex- 
tended cryostability model described by Eq.([T2]) and as- 
sume a power-law conductor with a voltage-current char- 
acteristic described by Eq.([T]). In other words we consider 
the helium temperature as variable and assume a smooth 
transition from the superconducting to the normal state. 
Additionally and opposite to other models [13, [ill we 
will consider that the index n is a function of the critical 
current and not a constant. The model equations with 
the explicit temperature dependence of Ic are 
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a system of two equation, one partial-differential and 
the other one pure algebraic. This system can be re- 
duced to a single partial-differential equation as follows. 
We use the first equation in ([T5|) to express the helium 
temperature as a function of conductor temperature 
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(16) 



and then we substitute it in the second equation. Cal- 
culating the derivative dThe/ dx is the only difficult part. 



First let us observe that dTcond/dx cannot be zero 
since for a power-law conductor G is never zero. Sec- 
ondly, it is clear physically that dTcond/dx should be pos- 
itive. The conductor temperature increases from some 
lower value at the inlet of the modeled section and reaches 
a maximum at its end due to the heat accumulated in the 
helium which is then convectively transported away along 
the conductor. Therefore a nontrivial solution of Eq. p^ 
is possible only if the following condition is satisfied 
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expressing G as Elq where Iq is the 
quench current, we arrive at the stability condition 
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(20) 

i.e the conductor with current dependent power-law 
index is stable against thermal perturbations if the local 
electrical field is smaller than Eq the take-off or quench 
electrical field. 

From experimental data we know that a — dn/dle > 
and /3 = dIe/dTeond < 0. The parameters a and /3 are 
slowly varying functions of temperature and field and 
can be considered in the first order of approximation as 
constants. This is true in particular for the linear depen- 
dence of n on le which is a good approximation for high 
temperature and/or fields i.e. for small critical currents. 
For constant n , Eq. ([20|) simplifies to 
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Expressing the power-law relation, Eq.([l] at the 
quench point: E = Eq at I = Iq 



Eq- Ee[j- 



one arrives finally at 



(22) 



Finally substituting Eq.(fT7|) in Eq. (|T5|) and using the 
first equation again we get finally 
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\nEeP, 

a relation obtained previously in [l^, [ll|. Equation 
(|20t is therefore the generalized expression for the quench 
condition in a conductor with variable power-law index 
n 



an equation for the conductor temperature Teond alone. 
This can be integrated numerically and the conductor 
temperature profile can be determined. This is however 
not our goal here. We are interested here much more 
in finding out the conditions under which a solution of 
Eg. dTSl) does exists. 



IV. RELATIONS TO THE EXPONENTIAL 

FORM OF THE VOLT-AMPERE 

CHARACTERISTIC 

An alternative form of the volt-ampere characteristic, 
used mainly in the Russian literature on this subject, is 



given by the following exponential form [ij 
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where Tj, and i?;, are a reference temperature and a ref- 
erence magnetic field. The reference point is chosen such 
that E = E,i,t Tcond = Tb,B = Bb and I = I^ (n, Bb). 
The constants Iq ,To and Bq are growth parameters of E 
by the current, temperature and magnetic field [14]. It is 
implicit assumed, but never explicitly stated, that these 
parameters are really constant i.e. they do not depend 
on temperature, current and field. The whole dependence 
on temperature, current and field is exactly what is seen 
in Eq. ([M|) explicitly and nothing else more. 

It the following we will establish the relation between 
this form and the power-law form of Eq.([T]) which for 
convenience is written here once again 



E {Tcond,B,I) — Ec 
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First let us observe a fundamental difference between 
the two functional forms. In the exponential form, the 
temperature and magnetic field dependence is explicit 
while in the power-law form it is implicit, through the 
temperature and field dependence of Ic- The connection 
between the two forms can be found by taking the first 
derivative of Ea.([M| with respect to Tcond at constant 
current and magnetic field. The result is 
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showing that indeed 1/Tq is the growth factor of an 
exponential grow. The same derivative calculated from 
Eq.dlSl) gives 
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(27) 



Now, as it is well known, the numerical difference 
between these two functional forms is vanishing small 
within 2 — 3 orders of magnitude and therefore we can 
equate the results from Ea. (P5)) and Eq. ipT)) with the final 
result 
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showing that the growth parameter is related to the 
slope of the critical current as a function of temper- 
ature and of n-index as a function of critical current. 
From Eq. ipS]) it is easy to see that a variable n makes 



To an explicit function of current and an implicit func- 
tion of temperature and field through the dependence on 
Ic- This contradicts the assumption on which expansion 
in Eq. ([M| is based i.e. that Tq should be a constant 
and not a function of anything. Therefore the expo- 
nential functional form is not compatible to the power- 
law with a variable index. For constant n, i.e. for 
dn/dic — and at a fixed temperature and field such 
that also did dTcond < is constant, we get an expres- 
sion which is independent of current 
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(29) 



a first relation between the parameters of the two func- 
tions. As can be seen the temperature growth parameter 
is connected to the n-index and the logarithmic deriva- 
tive of the critical current. A similar relation is obtained 
immediately for I = Ic from Eg.pS]) without any other 
supplementary conditions. 

For the magnetic field dependence we get using the 
same procedure, i.e. equating the derivatives dE/dB of 
Eq.dMI) and Eq.ipSj) 
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which is again current dependent through the logarith- 
mic term. For constant n or at I ^ Ic we get 
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Finally, for the current growth factor the following re- 
lation is obtained 
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It can be seen that in principle all growth parameter 
violate the assumption of complete separability of tem- 
perature, current and field dependence if the index of the 
power law depends on temperature and field. Therefore, 
in general, the exponential form is incompatible with 
a power-law with variable index. However, at I = Ic 
and/or at fixed temperature and field the two expres- 
sions are in agreement. 

The expressions for the growth parameters at I = Ic 
can be used to re-express the exponential form. Substi- 
tuting Eqs.(I2ni) ,(1311) and ^ in Eq.dMl we obtain 
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where h = Ic{Tb,Bt) and Ic = Ic{T,B). Observing 
now that the last three terms in the exponential repre- 
sent the Taylor series development of Ic {T, B) around the 
point (Tb, Bfy) we get an expression for the exponential 
form where the power-law index appears explicitly 



E — Ec exp 
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(35) 



which can be used as an alternative to the standard 
exponential form Eq. (1241) if the power-law index depends 
on temperature and field. Further, it can be shown that 
the new expression is equivalent to a power-law function 
to order 0{\/-n?). For this purpose we first invert Ea. ((35| 



where on the right-hand side we have developed the 
function in powers of \/n up to the second order. The 
logarithmic function is a slowly varying function and for 
large n, we see immediately that to order l/n^, Eq. (j36p 
and Eq. ([57|) are identical. 

In conclusion, the power-law form and the exponential 
form are identical to order 0(l/n^). For applications 
where the temperature and field dependence of the n- 
index is essential, the form in Eg.psp is recommended, 
whereas for currents close to critical current, the standard 
parametric form, Eg. ipS]) . can be used. 

The quench condition expressed with the help of the 
exponential form can be easily obtained from Eq. (|21[) and 
Eq.i 
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The inverted form of the power-law function is 
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a very simple and elegant relation between the local 
electric field at quench, the quench current and the tem- 
perature growth parameter. However, the applicability of 
this equation is limited to cases with constant or slowly 
varying n as discussed before. 
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